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Abstract 

Quantum operations provide a general description of the state changes al- 
lowed by quantum mechanics. The reversal of quantum operations is impor- 
tant for quantum error-correcting codes, teleportation, and reversing quantum 
measurements. We derive information-theoretic conditions and equivalent al- 
gebraic conditions that are necessary and sufficient for a general quantum 
operation to be reversible. We analyze the thermodynamic cost of error cor- 
rection and show that error correction can be regarded as a kind of "Maxwell 
demon," for which there is an entropy cost associated with information ob- 
tained from measurements performed during error correction. A prescription 
for thermodynamically efficient error correction is given. 
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I. INTRODUCTION 



Quantum operations arise naturally in the study of noisy quantum channels, quantum 
computation, quantum cryptography, quantum measurements, and quantum teleportation. 
In each of these applications it is of interest to learn when a quantum operation can be 
reversed. This paper gives a simple, physically meaningful set of necessary and sufficient 
conditions for determining when a general quantum operation can be reversed, thereby 
unifying and extending earlier work [1]-|3|]. The picture we develop applies equally well 
to the reversal of quantum measurements, as in the processes described by Mabuchi and 
Zoller H, and to the protection of quantum states against decoherence, as described in 
the literature on quantum error correction (see, for example, || for references). Finally, 
teleportation || can be understood as the reversal of a quantum operation, as was shown 
in ||, and the results obtained here are being applied in further work on characterizing 
schemes for teleportation. 

The paper is organized as follows. In Sec. we review the formalism of quantum 
operations and its application to the theory of generalized measurements, and we define 
the notion of a reversible operation. Section |TJ introduces information-theoretic measures 
associated with a quantum operation, which are used throughout the remainder of the paper. 
These measures, entanglement fidelity and entropy exchange, were introduced earlier by one 
of us 0] for the special case of deterministic operations; their definitions and properties are 
extended to general operations in Sec. fTJ. Section [IV A| states and proves an information- 
theoretic characterization of reversibility for general quantum operations, which is then given 
some simple applications in Sec. [IV B| . Section presents an alternative algebraic description 
of when a quantum operation can be reversed. Though many of the results in Sec. 
are already known, we provide new proofs, and the constructions used in these proofs are 
important in Sec. [VT] , where we give a thermodynamic analysis of error correction and show 
that schemes for performing perfect error correction can be done in a thermodynamically 
efficient way Section |V11] offers concluding remarks. 



II. QUANTUM OPERATIONS 

A. Definition and characterizations 

A simple example of a state change in quantum mechanics is the unitary evolution 
experienced by a closed quantum system. The final state of the system is related to the 
initial state by a unitary transformation U, 

p - £{p) = Upl/ï . (2.1) 

Unitary evolution is, however, not the most general type of state change possible in quantum 
mechanics. Other state changes, not describable by unitary transformations, arise when a 
quantum system is coupled to an environment or when a measurement is performed on the 
system. 

How does one describe the most general possible quantum-mechanical dynamics that 
takes input states to output states? The answer to this question is provided by the formalism 
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of "quantum operations." This formalism is described in detail by Kraus and is given a 
short, but quite informative review in an appendix to Q. In this formalism the input state 
is connected to the output state by the state change 

S(p) 

9 tr(£(p)) ' {2 - 2} 

which is determined by a quantum operation S. The quantum operation is a linear, trace- 
decreasing, completely positive map. Trace decreasing means that tr(£(p)) < 1 for all 
normalized density operators p. Complete posítivity means that in addition to preserving 
the positivity of density operators, the map preserves the positivity of all purifications 
of density operators. The trace in the denominator is included in order to maintain the 
normalization condition tr(p) = 1. 

The most general form for a completely positive map £ can be shown to be [[?].§ 

£(p) = £^Mj. (2.3) 

j 

The system operators Aj, which must satisfy 

£M-< J (2-4) 

j 

in order that 8 be trace decreasing, completely specify the quantum operation. We call 
Eq. ( |2.3p an operator-sum decomposition of the operation, and we refer to the operators Aj 
as decomposition operators. 

The operator-sum decomposition for a quantum operation is not unique, in that another 
set of decomposition operators {Aj} can give rise to the same operation. For example, the 
operation on a spin-| system defined by 

E,(p) = -Lp-L + ^-p^- (2.5) 

y/2 y/2 y/2 y/2 

can also be written in the form 

ei \ 1 + °z I + cr z I -a z I -a z 

£(P) = — P— + — P— ■ (2.6) 

Choi H has classified all sets of decomposition operators that give rise to the same 
operation. The result is that two sets of decomposition operators, {Ak} and {Bj}, give rise 
to the same quantum operation if and only if they are related linearly by a square unitary 
matrix u: 

Bj = Y J u jk A k . (2.7) 

k 

It is generally necessary to add some zero decomposition operators to the set with the smaller 
number of elements so that both sets have the same number of decomposition operators. We 
call a decomposition minimal if no decomposition into a smaller number of operators exists; 
a decomposition is minimal if and only if the operators in the decomposition are linearly 
independent. 
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We say that an operation 8 is pure if it can be written in terms of an operator-sum 
decomposition that contains only one operator; that is, there exists an operator A such that 



£{p) = ApA j . (2.8) 

The unitary transformation fl2.1|) is an example of a pure quantum operation. 

We say that an operation T> is deterministic or trace-preserving if tr(Z>(p)) = 1 whenever 
the input is a normalized density operator p. For a deterministic operation, the decomposi- 
tion operators Dj satisfy a completeness relation 

(2.9) 

j 

which implies that tr(P(p)) = 1. Notice that a pure deterministic operation must be a uni- 
tary transformation. Any deterministic quantum operation T> can be obtained by adjoining 
an ancilla system to the system of interest, allowing the system plus ancilla to interact 
unitarily, and then discarding the ancilla. Such a dynamics leads to a state change of the 
form 

p -> tr A (V(p ® o A )VÏ) = V(p) , (2.10) 

where tr A denotes tracing out the ancilla, a is the initial state of the ancilla, and V is the 
unitary operator for the joint dynamics of the system and ancilla. 
For a general quantum operation, 

ti(£{p)) = ti(p 1 £A]A J ) (2.11) 



is generally less than one, so S(p) must be renormalized, as in Eq. (|2T2|) , to produce an output 
density operator. A general quantum operation cannot be represented solely in terms of the 
joint unitary dynamics of the system and an ancilla, for that always leads to a deterministic 
operation, as in Eq. ( |2.1Ü| ). A general quantum operation can be obtained, however, if 



the joint dynamics is followed by a measurement on the ancilla; the quantum operation 
corresponds to particular measurement outcome described by an ancilla projection operator 
P A . The resulting state change, once the ancilla is discarded, is given by |7||| 

^ tvA{{i®P A )v{p®° A )v ] {i®P A )) £(p) f2 12) 

P ^ ti((I ®P A )V(p®a A )Vï{I ®P A )) tr(£(p)) ' { ' ' 

Notice that tr(£(p)) is the probability of the measurement result described by P A . A 
deterministic operation arises in the special case P A = I A . We note that for any quantum 
operation, there is a representat ion of the form ( |2.12 ) in which the initial ancilla state a A is 
a pure state. 

Suppose the measurement on the ancilla is described by a complete set of orthogonal 
projection operators P A , where the index i labels the measurement outcomes. Outcome i 
corresponds to a quantum operation 

£i(p) = tr A ((I ® P A )V{p ® <J A )V\I ® P A )) = tv A ((I ® P A )V{p ® a A )V r ) , (2.13) 
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which gives the unnormalized post-measurement state of the system, conditioned on outcome 
i. The probability for result i is 

Pi = tT(£i(p)) . (2.14) 

If one discards the measurement outcome, the output density operator is obtained by averag- 
ing over the outcomes, and the state change is given by a deterministic quantum operation: 

'^«sg&rÇ^- 1 ™- (2 ' 15) 

Thus a deterministic quantum operation can always be regarded as describing a measurement 
on the ancilla, whose result is discarded. 

In the case of a deterministic operation, the ancilla can be regarded as the system's 
environment; interaction with the environment gives rise to the nonunitary evolution that 
is described by the operation. For a general operation, the ancilla can also be regarded as 
an environment that can be observed and thus that has features of a measuring apparatus. 
Keeping these connotations in mind, we use the terms ancilla and environment interchange- 
ably in the remainder of the paper. 



B. Operations and generalized measurements 

The connection of quantum operations to quantum measurements is easy to explain. 
Standard textbook treatments describe quantum measurements in terms of a complete set 
of orthogonal projection operators for the system being measured. This formalism, however, 
does not describe many of the measurements that can be performed on a quantum system. 
The most general type of measurement that can be performed on a quantum system is 
known generalized measurement [f7| JX0|JXÏf1 . 



Generalized measurements can be understood within the framework of quantum oper- 
ations, because any generalized measurement can be performed by allowing the system to 
interact with an ancilla and then doing a Standard measurement described by orthogonal 
projection operators on the ancilla. Thus, as we can see from Eq. Q2.13|) , the most general 
type of quantum measurement is described by a set of quantum operations where the 
index i labels the possible measurement outcomes. The sum of the operations for the various 
outcomes is required to be a deterministic quantum operation, as in Eq. ( 2.15Q . 

Since we can give an operator-sum decomposition for each operation, 

£(p) = $>^4 , (2.16) 

j 

we can also say that the generalized measurement is completely described by the system 
operators Ay, which are labeled by two indices, i and j, and which satisfy the completeness 
relation 

E4^ = / ' ( 2 - 17 ) 
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If result i occurs, the unnormalized state of the system immediately after the measurement 
is given by 

£í(p) = J2 A vP A tj • ( 2 - 18 ) 

i 

The probability for result i to occur is 

p t = tr(£(p)) = tr^AljAj) ■ ( 2 - 19 ) 

This form makes the connection to the formalism of positive-operator-valued measures. The 
operators 

E^^A^Aj (2.20) 

3 

are elements of a decomposition of the unit operator into positive operators, as in the 
completeness relation ( |2.17|) . Such a decomposition of unity is called a positive-operator- 
valued measure (POVM). 

We say a measurement is pure if for each measurement result i, the corresponding quan- 
tum operation Ei is pure; that is, there exist operators Ai such that 

£i(p) = ApAl . (2.21) 

The probability that result i occurs is given by 

Pi = tr(pAtA) , (2.22) 

It can be shown that pure measurements correspond to extracting the maximum amount 
of information about the system from the state of the apparatus to which the system is 
coupled. 



C. Reversal of a quantum operation 

When we talk about reversing a quantum operation £, we generally do not mean that £ 
can be reversed for all input states, but rather only that £ can be reversed for all input density 
operators p whose support lies in a subspace M of the total state space L. In the case of a 
trace-preserving operation £, the subspace M is sometimes called a quantum error- correcting 
code or simply a code. It makes sense to talk about reversing £ on a subspace M only if 
£{p) 7^ for all p whose support lies in M, and we assume this condition henceforth. We 
say that a quantum operation £ is reversible on a subspace M if there exists a deterministic 
quantum operation 71, acting on the total state space L, such that for all p whose support 
lies in M, 

( m ) no£( P ) 
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Here 1Zo£ denotes the composition of 1Z with S, that is, TZo£{p) = 1Z(£(p)). We require the 
reversal operation 1Z to be deterministic because we want the reversal definitely to occur, 
not just to occur with some probability, conditional on some measurement result or ancilla 
state. 

In the problem of reversing deterministic quantum operations was considered. This 
case is of particular interest in situations where one is unable to obtain information about 
the environment. In contrast, B and M considered reversal of operations representing 
measurements, in which case information about the environment is available. 

We say that a measurement is reversible on a subspace M of the total state space L if for 
each measurement result i, the corresponding quantum operation is reversible. Outcomes 
that have zero probability on M are irrelevant, because they never occur, and thus they can 
be discarded. We could define measurements that are only sometimes reversible by requiring 
that only some of the measurement results have reversible quantum operations. Although we 
do not deal explicitly with such sometimes reversible measurements in this paper, the results 
obtained in Secs. [IV A| and |V A| , since they are derived for individual quantum operations, 
apply to sometimes reversible measurements. 



III. INFORMATION-THEORETIC MEASURES FOR QUANTUM OPERATIONS 

Schumacher [||] introduced entanglement fidelity and entropy exchange as useful informa- 
tion-theoretic measures for characterizing deterministic quantum operations. This section 
extends to general quantum operations the defmitions of entanglement fidelity and entropy 
exchange and generalizes the properties of those quantities obtained in JO and in [Q]. We 
begin by outlining the particular method for characterizing quantum operations that was 
used in JL] and [0 and that we use throughout the remainder of this paper. 

A. Method for characterizing quantum operations 

Suppose we have a quantum system, denoted henceforth by Q, and a quantum operation 
S that acts on states of Q. We denote the dimension of the Hilbert space of Q by D. It 
is convenient to introduce two mathematical artifices, a reference system R, whose Hilbert 
space has the same dimension, D, as the Hilbert space of Q, and an environment E, which 
has a Hilbert space of arbitrary dimension. 

The joint state of the system Q and the reference system R is chosen so as to purify the 
initial state of Q; that is, RQ is initially in a pure state p RQ = \ïf/ R Q)(ïff R Q\ satisfying 

tr R (\y RQ )(V RQ \) = p Q , (3.1) 

where p® is the initial state of system Q. To reduce the clutter in the notation, we drop 
the Q superscript when it is clear that we are dealing with the primary quantum system Q. 
The initial state of the environment E is assumed to be a pure state p E = |e)(e|, which is 
uncorrelated with the system RQ. Thus the initial state of the overall system is also pure: 
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The joint system QE is subjected to a two-part dynamics consisting of a unitary oper- 
ation, U® E , followed by a projection onto the environment alone, described by a projector 
P E . The reference system R has no internal dynamics and does not interact with Q or E. 
This two-part dynamics leaves the overall state pure. As we mentioned above, it is always 
possible to find \e), U QE , and P E such that 

S(pQ) = tr E (P E U QE (p Q ® p E )U QE] P E ) = tr RE (P E U QE (p RQ <g> p E )U QE] P E ) . (3.3) 

We denote the normalized states of the different systems R, Q, and E after this evolution 
by primes. Of special interest is the joint state of RQ after the dynamics, which is given by 

RQI _ tr E (P E U QE (p RQ g p E )U QE] P E ) (l R ®S)(p R Q) _ (l R ®S)(p RQ ) 
9 te{P E UQ E {p R ® ® P E )UQ E ^P E ) ~ ti((l R ® £)(p RQ )) ~ tr(£(jfl)) ' 

where X R is the identity operation for the reference system. Using the operator-sum decom- 
position ( |2.3| ) of £, we can write p R ® in terms of an ensemble of unnormalized pure states, 
that is, as a sum of terms each of which is proportional to a one-dimensional projection 
operator: 

rq' _^ (I R ® A)) 

The state of the reference system after the dynamics is given by 

tx{P E UQ E {p R Q ®p E )U^P E ) QK ' 

We emphasize that p R ' is generally not the same as p R , because of the presence of the 
environment projector P E in Eq. (|3.6| ). This is in contrast to the case of a trace-preserving 
operation S, for which the environment projector is absent, that is, P E = I E in Eq. ( |3.6| ) 
and, hence, for which p R = p R . 



B. Entanglement fidelity and entropy exchange 

Following Schumacher 0,0] , we define the entanglement fidelity to be the fidelity with 
which the joint state of RQ is preserved by the dynamics: 

f . fag)B <^^|^ = <*!g!gWm . (3.7, 

Using the form (|3.5| ) of p R ®' and noting that 

(q R Q\(I R ® Aj)\^ RQ ) = tr(\^ RQ )(^ RQ \(I R (8) Aj)) = tr (p Q A 3 ) , (3.8) 

we can put the entanglement fidelity in a form that, as implied by our notation F e (p,S), 
manifestly depends only on the initial state of Q and the operation that is applied to Q, 
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(3.9) 



The entropy exchange is defined to be 

S e (p,£) = S(p RQ ') = S(p E ') , (3.10) 

where S(p) = — tr(plogp) denotes the von Neumann entropy of the density operator p 
and where the latter equality follows from the fact that the overall state after the two- 
part dynamics is pure. This generalizes the definition of entropy exchange in ffl to general 
quantum operations £ . 

The entropy exchange obeys several inequalities that follow from the subadditivity of von 
Neumann entropy |[2|-l4}| and the purity of the overall state of RQE after the dynamics: 



S e = S(p R V) < S(p R ') + S(pQ') , (3.11) 
S(p Q ') = S(p RE ') < S(p R ') + S(p E ') = S(p R ') + S e , (3.12) 
S(p R ') = S(p^ E ') < S(pQ') + S(p E ') = S(pV) + S e . (3.13) 

Each inequality here is an expression of subadditivity, with equality holding if and only if 
the joint density operator on the left factors into a product of the two density operators on 
the right (for example, equality holds in Eq. Q3.11| ) if and only if p R ® = p R ® p® ). The last 
two of the above inequalities can be combined into a single inequality, 

S e (p,£) = S(p R V) > \S(ffl) - S(p R ')\ , (3.14) 



sometimes known as the Araki-Lieb inequality 

If £ is trace preserving, then as noted below Eq. (|3.6Q , the state of the reference system 
does not change under the dynamics, that is, p R = p R ; moreover, since the initial state of 
RQ is pure, we always have that S(p R ) = S(p Q ). Thus, for a trace-preserving operation, 
we can use S(p R ) = S(p®) to eliminate the reference system from the above inequalities, 
leaving the inequalities obtained by Schumacher For a general quantum operation, it is 
not true that S(p R ) = S(p Q ), and the inequalities must be left in the form given above. 

The von Neumann entropy of p R ® is the same as the entropy of the matrix of inner 
products formed from the unnormalized pure states that contribute to the ensemble for p R ® 
in Eq. ( |3.5| ). Explicitly, in terms of a positive, unit-trace matrix W, whose components are 

(^\(I R ®Al)(I R ®A,)\^) 

w > k ~ WÏp^ï) ' ( } 

the entropy exchange is given by 

S e (p,£) = S(W) = -tr(WlogW) . (3.16) 
The components of W can be simplified to the form 

_ tr(\fy RQ )(fy RQ \(I R g A\)(I R g Aj)^ _ tr(A lP A{) 
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Notice that the diagonal elements of W, 



* = w - = i£é • (3 ' 18) 



are the probabilities with which the pure states in Eq. Q3.5|) contribute to the ensemble. 

Relative to a particular density operator p, there is a "canonical decomposition" of the 
quantum operation £. Suppose the matrix W arises from p and a particular operator-sum 
decomposition of £ in terms of decomposition operators Aj, as in Eq. ( J3.17D , and suppose 
we diagonalize W with a unitary matrix u, 

S U jl W lrnU*km = X j 5 jk , (3-19) 
l,m 

where the nonnegative real numbers Xj are the eigenvalues of W. Now define new operators 

A j = Y i 'u·3kA k . (3.20) 

k 

These operators being a unitary remixing of the original decomposition operators, they give 
another operator-sum decomposition for £, with an associated matrix 

W jk = tr t ^"(ffi = E UiiWimU* km = XjS jk . (3.21) 



We say that a decomposition Aj satisfying Eq. ( p.21|) is a canonical decomposition of 



£ with respect to p. The canonical decomposition is unique up to degenerades in the 
eigenvalues Xj and up to (trivial) phase changes in the canonical decomposition operators 
Aj. The entropy exchange can be written as 

S e = S(W) = S(W) = - y £X j logX j =H(X), (3.22) 

j 

where H(X ) is the Shannon information constructed from the probability distribution À. It 
can be shown that 

Se<H(?) = - y Eq j logq j , (3.23) 

j 

equality holding only for a canonical decomposition. 

Notice that for a pure quantum operation, the decomposition that contains only a single 



decomposition operator, as in Eq. Q2.8D , is the canonical decomposition with respect to any 
density operator p. The canonical W matrix is the one-dimensional unit matrix, and hence 
the entropy exchange S e is zero. 

Consider the canonical decomposition of an operation £ with respect to the unit density 
operator I/D, where D is the dimension of the system Hilbert space. Such a canonical 
decomposition, whose decomposition operators satisfy 

tr(ÀjÀl) = Xjtr(£(I))6 jk , (3.24) 
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is a minimal decomposition of £. All minimal decompositions can be obtained from this one 
by unitary remixings that leave the number of operators in the decomposition unchanged. 
The entanglement fidelity and the entropy exchange obey the quantum Fano inequality, 

S e (p, £) < h(F e (p, £)) + (1 - F e (p, £)) \og(D 2 - 1) , (3.25) 

where h(p) = — plogp — (1 — p) log(l — p) and D is the dimension of the system Hilbert 
space. The quantum Fano inequality was first derived by Schumacher [||] for trace-preserving 
operations, but Schumacher's proof goes through unchanged for general quantum operations. 



C. Data-processing inequality 

In the following we often consider a two-step operation that consists of two successive 
operations. The situation of interest here is that of reversing an operation, as discussed in 
Sec. [II Q : the first of the two operations is the operation to be reversed, and the second 
is a reversal operation, which is necessarily deterministic. Schumacher and Nielsen de- 
rived an important inequality, called the data-processing inequality, for the case in which 
both operations in a two-step operation are trace preserving. Here we show that the data- 
processing inequality remains vàlid for any two-step operation in which the first operation 
£ is arbitrary, but the second operation T> is deterministic. 

In this situation we use E to denote the environment used in the first step and A 
to denote the ancilla or environment used in the second step. We use double primes to 
distinguish the normalized states of the various systems after the second step. We can 
draw several conclusions about the von Neumann entropies of various states in this scenario. 
In particular, since the overall state after both steps is pure, we have that S(p REA ) = 
S(p Q ) and that S(p EA ) = S(p RQ ) = S e (p, Vo£) is the entropy exchange for the two-step 
operation. Moreover, since the state of RQE is pure after the first step and since R and E 
do not participate in the second step, we have that S(p RE ) = S(p RE ) = S(p Q ) and that 
S(p E ) = S(p E ) = S(p R ® ) = S e (p,£) is the entropy exchange in the first step. 



The strong subadditivity property of von Neumann entropy [13—15] constrains the en- 
tropies after the two-step dynamics: 

S(p REA ") + S(p E ") < S(p RE ") + S(p EA ") . (3.26) 

Substituting the entropy relations just derived and re-arranging yields the data-processing 
inequality, 

S(p Q ') - S e (p, £) > S(p Q ") - S e (p, Vo£). (3.27) 

The left-hand side of the data-processing inequality is constrained by Eq. (|3.12| ), leading to 
the double inequality 

S(p R ') > S(pV) - S e (p,£) > S{pQ") - S e (p,Vo£) . (3.28) 

If £ is trace preserving, then S(p R ') = S(p®) and this double inequality reduces to the form 
found by Schumacher and Nielsen 0]. 
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IV. INFORMATION-THEORETIC CHARACTERIZATION OF REVERSIBLE 

QUANTUM OPERATIONS 



A. General information-theoretic characterization 



In this section we demonstrate that a general quantum operation £ is reversible on a 
subspace M of the total state space L if and only if the following two conditions are satisfied: 

Condition 1: tr(£(p)) = p 2 for all p whose support is confined to M, (4-1) 

where /i is a real constant satisfying < p < 1; 

Condition 2: 

S(p) = SÍ — ) — SJp,£) for any one p whose support is the entirety of M, 
\tï{S{p))J 

(4.2) 

(and then for all p whose support is confined to M). 
Condition 1 is equivalent to 

P M EP M = P 2 Pm , (4.3) 

where 

E = J2A]Ai (4.4) 

j 

is the POVM element corresponding to £ and Pm is the projector onto M. Condition 1 
has the appealing intuitive interpretation that if we view £ as a dynamics for the system, 
conditional on some measurement result or post-interaction environment state, knowledge 
of that result or state gives no information about the initial system state p. Condition 2, 
though less intuitive, states essentially that for initial states whose support lies in M, no 
quantum information is lost in the dynamics described by £. 

We begin by proving necessity. Suppose that £ is reversible on M. Then it was shown 
in || that Condition 1 follows. The reason is that reversibility implies that 1Z o £ (p) = 
tr(£(p))p for all p whose support is confined to M; 1Z o £ being linear, the only way this 
equation can be satisfied is if tr(£(p)) has a constant value p? > 0. Let £u be the restriction 
of £ to M, that is, 

£ M (p) = ]T AjPmpPmA] . (4.5) 
j 

Notice that £m(p) = £{p) if p has support lying wholly in M. Let M be the subspace 
that is the orthocomplement of M and Pj^ be the projector onto M. Now introduce a new 
quantum operation J 7 , whose action on any p is given by 

HP) - ^# + PJIPPJI ■ (4-6) 
p 
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The reason for introducing T is that it is a deterministic operation with the property that 



Tip) = £ M = JÏÍEL (4 7) 

for states p whose support lies wholly in M. Thus £ is reversible on M if and only if T is 
reversible on M. Since T is deterministic, however, the necessary and sufficient condition 
for its reversibility is the condition already obtained by Schumacher and Nielsen 0: T is 
reversible on M if and only if 

S(p) = S(f(p))-S e (p,f) (4.8) 

for any one p whose support is the entirety of M (and then for all p whose support lies in 
M) . We complete the proof of necessity by noting that for states p whose support is confmed 
to M, the W matrices of £ and T are the same, which implies that 

S e (p,£) = S e (p,F). (4.9) 



Substituting Eqs. fl4.9|) and (|4.7| ) into Eq. ( |4.8j ) yields the second condition ( |4.2| ) 



The sufficiency of Conditions 1 and 2 is proved in an obviously similar way, but one point 
should be stressed. For T to be reversible on M, it is sufficient that Eq. ( |4.8| ) hold for any 
one p whose support is the entirety of M. Thus for £ to be reversible, it is sufficient that 



the second condition ( |4.2|) hold for any one such p. 

Before going on, one further point deserves mention. If the initial state of Q is the unit 
density operator in the subspace M — that is, p = Pm/oI, where d < D is the dimension 
of M — then we can dispense with Condition 1. What we are claiming is the following 
equivalence: £ is reversible on M if and only if 

\ogd = S{P M /d) = s if ( ^f d)) ) - S e (P M /d,£) . (4.10) 

The necessity of Eq. ([4.101) has already been shown. We now demonstrate sufficiency by 
showing that Eq. ( 4.10| ) implies the first condition ( |4.1| ). 



For this purpose, notice that when p Q = Pul d, the initial pure state of RQ is an 
entangled state of the form 

I^H-^Ebd}®!^)- (4.11) 



m=l 



Here the kets are an orthonormal basis for the <i-dimensional subspace M, and the 
kets \Xm) are a se t °f d orthonormal vectors for R. Substituting this entangled state into 
Eq. ( p.5| ) yields a new expression for the joint state of RQ after the dynamics: 

RQ' = V- \ y R\(yR\ (g) 1 y M^m)^^ 
P L·lXmHXnl®^ tm p M , d)) ■ ^ 

The state of the reference system after the dynamics now assumes the form 
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n R' — f r í n RQ' \ — X^ ^ V/Vjjr_m / \R\/_,R\ (A-l-i\ 

p - tlQÍP ï-^dKMPM/d)) lxJ(xJ ' (413) 

where i? is the POVM element associated with £ [cf. Eq. (fO])]. Since the support of 
p R lies within the <i-dimensional subspace spanned by the vectors we know that 

S{p R ) < logd. Condition (|4.10| ), when combined with the left inequality in Eq. fl3.28| ), 



implies that S(p R ') = logd. This means that the matrix elements of p R ' in Eq. ( 4.13 ) must 
be those of the unit density operator Pm/ú on M. Hence we conclude that 

P M EP M = ti(S(P M /d))P M , (4.14) 

which as already noted in Eq. ( |4.3| ), is equivalent to Condition 1, with p? = ít{£(Pm/ d)). 



B. Applications of information-theoretic conditions for reversibility 

A number of useful results follow from the information-theoretic characterization of re- 
versible quantum operations found in the preceding subsection. Among these is a general 
characterization of teleportation schemes, which has been discussed in || and is the subject 
of continuing work. This subsection describes several simpler, but still useful applications 
of the information-theoretic characterization. 

We first show that an operation is reversible by a unitary operation if and only if it acts 
like a múltiple of a unitary operation when restricted to the reversal subspace. 

Theorem. A quantum operation 8, with decomposition operators Aj, is reversible by a 
unitary operator [/ona subspace M if and only if there exist complex constants Cj such 
that 

A 3 P M = c 3 U^P M . (4.15) 

The constants satisfy 

EM 2 = A ( 4 - 16 ) 

3 

where p? is the constant value of tr(£(p)) on M. 



Proof. The sufficiency of the condition ( 4.15|) is obvious. The proof of necessity is to 



notice that for all p, not just those whose support is confined to M, we have 

U^Mrf = P MP P M , (4.17) 

where Sm is the restriction of p to M. Rewritten as 

Mp) = E AjPupPuA] = p 2 U^P M pP M U , (4.18) 

j 

this shows that Sm is a pure operation whose canonical decomposition contains the single op- 
erator pWPm- By the result fl2.7|) that relates operator-sum decompositions, the conclusion 
follows. This completes the proof. 



14 



A second theorem shows that an operation that is reversible for all initial states acts like 
a múltiple of a unitary operation. 

Theorem. A quantum operation £ that is reversible on the entire state space of the 
system is a positive múltiple of a unitary operation; that is, 

S(p)= l· L 2 UpUï (4.19) 

for some constant p satisfying < p < 1. 

Proof. A simple proof can be obtained by examining the reversal operation constructed 
in and verifying that it is unitary. The result follows from this and the fact that p 2 = 
tr(£(p)) is a constant. 

We present here, however, a purely information-theoretic proof that does not require the 
explicit construction of a reversal operation. A quantum operation that is reversible on the 
entire D-dimensional state space of Q satisfies 



\ogD = S(I/D) = s(W^\ - S e (I/D,£) . 



(4.20) 



Since the entropy is maximized by I/D and the entropy exchange is nonnegative, we see 
immediately that S (£ (I / D) / fi 2 ) = logD and S e (I/D,£) = 0. The first of these conclusions 
means that 

£(I/D) = p 2 (I/D). (4.21) 

The second means that the W matrix for initial state I/D is of rank one; thus there exists 
a unitary matrix u such that 

Y, u 3i w i™ u % m = Sji6 kl . (4.22) 

l,m 

Defming a canonical decomposition of £ as in Eq. (|3.20| ), the canonical W matrix becomes 

W jk = = 5jl 8 hl . (4.23) 

It follows that only the first operator in the canonical decomposition, A\, is nonzero. Com- 
bining this result with Eq. ( (4.21| ) yields 

£{I/D) = À X À[/D = p 2 (I/D) , (4.24) 

which implies that A\ = // U. This completes the proof. 

This second result has a useful application to teleportation. Recali the bàsic set-up for 
teleportation |J. Alice possesses an unknown quantum state p, which she wishes to teleport 
to Bob. Alice also sends Bob some classical information, which we represent by i. It was 
shown in |§ that the state of Bob's system, conditioned on the information i, is related to 
the state of Alice's system by a quantum operation £ 4 , 

tr(íi(p)) 
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If Bob wishes to achieve teleportation then he must be able to reverse the operation £ What 
the above result shows is that Bob can use a unitary operation to do the reversal, since his 
reversal must work over the entire space of initial states p. No generality is introduced by 
allowing Bob to use nonunitary reversal operations, that is, by allowing Bob to employ an 
ancilla to assist in teleportation. Considering only unitary reversals, as was done in [[J, is 
thus sumcient for the study of teleportation. 

We can also compare the results obtained in this paper to earlier characterizations of 
reversible pure and deterministic quantum operations, obtained in || and || (where pure 
operations are called ideal operations), and show that these earlier results are special cases 
of the general characterization embodied in Conditions 1 and 2. 

For a deterministic quantum operation £, it is certainly true that tr(£(p)) = 1 is a 
constant, so Condition 1 is automàtic. Thus reversibility for a deterministic operation 
is equivalent to Condition 2 alone, that is, to S(p) = S(£(p)) — S e (p,£), which is the 
reversibility condition obtained in 0. 

For a pure quantum operation £, Condition 1 is equivalent to P M A^ AP M = p 2 Pm- This 
implies, using the polar-decomposition property of operators, that APm = \/ PmA^APm = 
pWPm for some unitary operator U. Hence £ can be reversed by U. Thus reversibility for 
a pure operation is equivalent to Condition 1 alone, that is, to tr(£(p)) = fi 2 , which is the 
reversibility condition obtained in ||. For a pure operation, Condition 2 can be dispensed 
with because it follows from Condition 1: for any pure operation the entropy exchange is 
zero, and Condition 1 implies that £ acts like a múltiple of a unitary on M, which means 
that S(p) = S(£(p)/p 2 ). 

V. ALGEBRAIC CHARACTERIZATION OF REVERSIBLE OPERATIONS 

Up to this point we have taken an information-theoretic approach to the reversal of 
quantum operations. In this section we switch to an algebraic point of view. The algebraic 
results obtained in this section are particularly powerful when used in combination with the 
information-theoretic viewpoint, as we illustrate in the next section on the thermodynamic 
cost of error correction. We begin with the theorem that establishes algebraic conditions for 
reversibility. 

A. Reversibility theorem 

Theorem. A quantum operation £ , with decomposition operators Aj, is reversible on M 
if and only if there exists a positive matrix m such that 

P M A\A 3 P M = m jk P M . (5.1) 

The trace of m, 

J2m jj = p 2 , (5.2) 

j 

is the constant value of tr(£(p)) on M. (Under a unitary remixing of the decomposition op- 
erators, the matrix m undergoes a unitary transformation, which leaves the trace invariant.) 



16 



This result was proved by Knill and Laflamme [16] and by Bennett et al. \ï%. We give 



a different proof, particular ly of the sufíiciency of condition (5J.). The construction used in 
our proof is crucial to our subsequent analysis of the thermodynamics of error correction. 

Proof. We deal first with the necessity of condition ( |5.1|) and notice that for all density 
operators, not just those whose support is confined to M, we have 



K o S M (p) = E RiAjPmpPmAIRJ = p 2 P M pP M 



(5.3) 



1,3 



where the operators Rk make up an operator-sum decomposition for the reversal operation 
TZ. Equation ( |5.3|) means that TZ o £ M is a pure operation, whose canonical decomposition 
consists of the single operator pPu- By the result ( |2.7| ) that relates operator-sum decom- 
positions, we can conclude that there exist constants Cji such that 



RiAjP\i 



CjlPu ■ 



The constants satisfy 



El |2 2 

3,1 



(5.4) 



(5.5) 



Using the trace-preserving property of the reversal operation TZ, we can write 

PmAÍAjPm = Ç PmAÍrIr^Pm = Cjl c* k ^ P M = m jk P M , (5.6) 

where the matrix m = cé is manifestly positive. 

We now demonstrate that condition (|5.1| ) is sufficient for reversibility. Let «bea unitary 
matrix that diagonalizes m, that is, 



Ujimi n u* kn = djSjk 

l.n 



(5.7) 



where the nonnegative real numbers dj are the nonnegative eigenvalues of m. Relative to a 
new decomposition of S , defined by 



Aj = UjkAk 



(5.8) 



condition (|5.1|) becomes 



PmA\AjPm — djójkPhí ■ 



(5.9) 



The diagonal (j = k) elements of Eq. ( |5.9j ) imply, by the polar-decomposition property, that 
there exist unitary operators Uj such that 



AjP M = U^PmA)A ó P m 



'djUjPj 



(5.10) 



Notice that if dj = 0, the corresponding decomposition operator Aj is irrelevant to the 
operation of £ within M, although such an Aj is generally important to the action of £ on 
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density operators whose support is not confined to M. When there are such decomposition 
operators, i.e., when dj = for some j, the subspace M is called a degenerate code; we 
discuss the meaning and significance of degenerate codes in Sec. |V B . For dj ^ we let Mj 
be the subspace that M is mapped to by Uj, and we let Pj = l/jP^Uj be the projector onto 
Mj. The off-diagonal elements of Eq. fl5.9| ) imply that these subspaces are orthogonal, that 
is, 

P k P, = 5 jk Pj . (5.11) 

The action of £ on any density operator whose support is confined to M takes the 
following form in terms of the new decomposition: 

£(p) = Ç À 3 P MP P M À) = J2 d,jU 3 P MP P M U] = ]T djPjUjpUjPj . (5.12) 
i j j 

It is easy now to construct an operation that reverses £. Let N be the subspace that is the 
direct sum of the orthogonal subspaces M,-, and let N be the orthocomplement of N. The 
projector onto N is given by 

P W =I-P N = I- £ P 3 . (5.13) 

{jK¥o} 

Now we define the action of a putative reversal operation by 

K{ P )= £ UjPjpPjUj+PjïpPü . (5.14) 

{Mm 

This reversal operation is trace preserving, as required, since 

E l',l)'V) ■ ''-''x E P 3 + Pn=I, (5-15) 

{j\d^0} {j\d,m 

and simple àlgebra shows that for all p whose support is confined to M, 

Ko£(p)= l· i 2 p, (5.16) 

where 

^ = J2d j =tr(£(p)). (5.17) 

i 

Thus 1Z is indeed a reversal operation for £ on the subspace M. This completes the proof. 



B. Discussion of algebraic conditions for reversibility 

The proof has a compelling physical interpretation. It shows that an operation that is 
reversible on M has an operator-sum decomposition in which the decomposition operators 
A 3 map M unitarily to orthogonal subspaces Mj. The operation on the entire space is gen- 
erally not representable as an ensemble of unitary operations, but as far as its action on the 
reversible subspace M is concerned, the operation can be represented by unitary operators 
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Uj, which are applied randomly with probabilities Aj = dj/p 2 and which, moreover, take 
M to orthogonal subspaces Mj. Reversal can be effected by first measuring in which of 
the orthogonal subspaces the state lies after the operation, thus determining which unitary 
operator Uj occurred, and then applying the corresponding inverse unitary operator Uj to 
restore the initial state. We stress that one can always effect reversal in this way, by using a 
measurement — indeed, a pure, projection-valued measurement — followed by a unitary con- 
ditioned on the result of the measurement. It is equally important to emphasize, however, 
that the deterministic reversal operation can also be constructed without measurements, by 
using an ancilla as described in Sec. [IT]. The two methods of reversal lead, of course, to the 
same reversal operation. 

The operator-sum decomposition Aj used in the above proof is obviously quite special. It 
is a canonical decomposition for £ relative to the initial state Pm / d, where d is the dimension 
of M, as can be seen directly by taking the trace of Eq. (|5.9|) . More interesting is that the 
operators Aj are a canonical decomposition for £ relative to any initial density operator p 
whose support lies in M; that is, the W matrix is diagonal, 



w tr (-v4) ti(AjP MP p M À\) dj 

LL LL LL 



with eigenvalues Xj. It follows that 

S e (p,£) = S(W) = H(X) . (5.19) 

Moreover, for any p whose support lies in M, Eq. ( |5.12| ) shows that the density operator 
after application of £ is given by 

f/ = e( P )/ f i 2 = y E*jPi, (5-20) 

where 

Pj = UjpU) (5.21) 

is the unitary image of p in Mj. It follows that S(pj) = S(p) and since the density operators 
Pj are orthogonal, that 

S(£(p)/p 2 ) = H(X) + S(p) = S(p) + S e (p, £) . (5.22) 



This is an explicit demonstration of condition ( |4.2|) . 

The existence of the canonical decomposition Aj of Eq. ( |5.8| ) clarifies the notion of a 
degenerate code. A common way of defining degeneracy is to say that a code is degenerate 
if any of the off-diagonal elements of the matrix m in Eq. ( |5.1| ) are nonzero. This defmition is 
flawed, however, because it is not invariant under changes in the operator-sum decomposition 
of £. The off-diagonal elements of m can always be made to vanish by transforming to the 
canonical decomposition Aj. 

Loosely speaking, what degeneracy is supposed to capture is the idea that some of the 
"errors" produced by £ are irrelevant within the code subspace M. This idea must be 
translated into a mathematical form that is independent of the operator-sum decomposition 
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of £. One way of doing so was introduced by Gottesman suppose that the operators A 



l 3 



constitute a minimal (and thus linearly independent) decomposition of E\ if the restricted 
operators AjP M , which form a decomposition of the restricted operation £ Ml are linearly 
dependent, Gottesman calls the code degenerate. The reason for this defmition is that if the 
operators AjPu are linearly dependent, then transformation to a minimal decomposition 
of Em reduces the number of decomposition operators, i.e., reduces the number of "errors." 
The canonical decomposition provides just such a minimal decomposition of Em-, that is, the 
decomposition consisting of the restricted operators AjPm- The reduction in the number of 
errors shows up in that the operators Aj that have dj = are irrelevant to the operation 
of E within M. Thus we arrive at the equivalent definition of degeneracy introduced in the 
above proof: a code is degenerate if one or more of the eigenvalues dj vanishes. 

This discussion leads to a manifestly invariant way of defining degeneracy in terms of 
operator subspaces: a code is degenerate if the operator subspace spanned by the decom- 
position operators of £ has higher dimension than the operator subspace spanned by the 
decomposition operators of Em- Moreover, it is now clear why degeneracy is considered a 
possible means of beating the "quantum Hamming bound." That bound is derived from 
counting the number of possible linearly independent errors, not restricted to the code sub- 
space, and assuming that error correction requires for each error an orthogonal subspace the 
same size as the reversible subspace M. 

We turn now to properties of the reversal operation. In Eq. Q5.14 ) we introduced a 



particular reversal operation 7Z, which is defined in terms of an operator-sum decomposition 
that consists of the operators 



and the operator 



Rj = ü}Pj = P M U] for j such that dj ^ 0, (5.23) 



Rn=Pn- (5-24) 



The important part of TZ is its restriction to the subspace N; the action of the restricted 
operation TZn is defined by 

n N (p) = Y,Rj p NpPNR] + RNPNpP N RN =Y,RjpR]= E UjPjpPjUj- (5.25) 

The first question we address is the extent to which the reversal operation is unique. 
For that purpose, consider another operation T, with decomposition operators ï], which 
reverses E on M. The action of T on an arbitrary density operator can be written as 

T(p) = £ TtPxpPxT? + £ TtPüpPjjT} + J2 T^pPjï + P„pP N )T} . (5.26) 
i i i 

The first and second terms in this expression are the restrictions of T to the orthogonal 
subspaces iV and N, respectively, and the third term is an additional contribution that can 
arise when p is not block-diagonal with respect to iV and N. The second and third terms 
are unaffected by the requirement that T be a reversal operation; the only restrictions on 
the second and third terms come from the requirement that T be trace preserving. The first 
term defines the action of the restriction of T to the subspace N, that is, 
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T N (p)=J2TiP N pP N T^ , (5.27) 
i 

The restricted operation is the important part of T for reversal. What we show now is 
that T N is the same operation as TZjy. 

We proceed by noting that the decomposition operators Ti must satisfy Eq. ( |5.4| ) for any 
decomposition of £ and, in particular, must satisfy it when the decomposition operators for 
£ are chosen to be the canonical decomposition Aj; that is, we must have 

JljTiUjPM = TiÀ^Pm = cjiPm , (5.28) 
where the constants Cji determine the diagonalized m matrix of Eq. ( |5.7|) , 

Y, d jicli = d j5jk ■ (5.29) 
i 

We now discard the vàlues of the index j for which dj = 0; this eliminates rows of zeroes 
from the matrix c. For the remaining vàlues of j, we have that 

TíPjUj = WjPm = ^LP M = v tj P M , (5.30) 
The columns of the matrix v are orthonormal, that is, 

S v ik v U = -/= £ è H è li = 5 rt (5-31) 
i ydjdk i 



(this means, in particular, that the number of rows of v is not smaller than the number 
of columns), and thus by adding columns, v can be extended to be a unitary matrix. By 
moving the unitary operator Uj in Eq. (|5.30|) to the other side of the expression, we obtain 



which implies that 



TiPj = v li P M ü) = v lj U}P ú = v^Rj , (5.32) 



T t P N = TiP 3 = £ VyRj . (5.33) 



Since the decomposition operators TiP N are related to the decomposition operators Rj by a 
unitary matrix, we can conclude, as promised, that 7/v and IZ^ are the same operation. 

The upshot is that the part of a reversal operation that actually effects the reversal — that 
is, the restriction of the reversal operation to the subspace N — is uniquely determined. In 
what follows this permits us to make general statements about all reversal operations. 

The decomposition used to define TZ in Eqs. ( |5.23|) and (|5.24|) is special. For any p whose 



support is confmed to M, this decomposition is a canonical decomposition for 1Z relative 
to the output state p' = £(p)/p 2 of Eq. (|5.20|) . This fact is crucial to our later analysis of 



the thermodynamic emciency of error correction. To prové it, notice that for any p whose 
support is confmed to M, we have Pjp'Pk = XjSjkPj and Pjqp' = p'Pjj = 0. It follows that 
the W matrix is diagonal, that is, 
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W jk = tr{Rjp'R\) = tr(UjP jP 'P k U k ) = \,5 jk (5.34) 



3 

and 

Wjn = W WNj = W W ü = . (5.35) 

The canonical W matrices for £ and 7Z being the same, the entropy exchange in the 
reversal operation is the same as the entropy exchange in £: 

S e (p',n)=H(\) = S e (p,S) . (5.36) 

In addition, since TZ(p') = p, Eq. ( |5.22j ) can be recast as 

s e (p',n) = s(p')-s(n( P ')). (5.37) 

This result, that the entropy exchange in the reversal operation is equal to the entropy 
reduction, is important for our discussion of the thermodynamics of error correction in 
Sec. [V|. Since the entropy exchange (|5.36|) is determined by the restriction of TZ to N, 



Eqs. ( |5.36| ) and ( |5.37| ) hold for operation that reverses £ on M. We stress, however, that 
Eq. (|5.37f) does not hold generally for trace-preserving operations; rather, as Eq. (|3.14|) 



shows, all that one can say for a general trace-preserving operation is that the entropy 
reduction does not exceed the entropy exchange. 

We can also make some powerful observations about the reversibility of entire classes of 
operations. Knill and LaFlamme |1| showed that if an operation £ , with decomposition 
operators A k , is reversible on M, then any operation J 7 , whose decomposition operators Bj 
are linear combinations of the A k , is also reversible on M. This can be seen immediately 
from Eq. Q5.1|) : if 

B j = Y / b 3k A k (5.38) 

k 

(since bj k is not assumed to be a unitary matrix, this is not just a unitary remixing, which 
would yield another decomposition of £ instead of a new operation), then 

P M B\B 3 P M = £ ÒI^PmAIAPm = ( £ b jimim bi m j P M = n ]k P M , (5.39) 

l,m \l,m J 

where the matrix n = bmtf is manifestly positive. 

This result can be stated compactly in the language of operator subspaces: the re- 
versibility of an operation £ on M implies the reversibility on M of any operation whose 
decomposition operators span a subspace of the span of the decomposition operators of £. 
We stress that the decomposition operators A k can be written as a linear combination of 
the decomposition operators Bj only if the operators Bj span the entire operator subspace 
spanned by operators Aj. 

We can go further to show that any operation that reverses £ is also a reversal operation 
for T . To do so, notice first that the decomposition operators Bj can be written as a linear 
combination of the canonical decomposition operators A k of Eq. ( |5.10| ): 
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k 

Thus the action of T on any density operator p whose support is confined to M can be 
written as 



T(p) = J2 BjP M pP M B) = hkb^ÀkPMpPMÀj = J2 ^J d kdi b úk %P k U k pU\P l , (5.41) 

3 j,k,l j,k,l 

The constant value of tr (jF(p)) is given by 

v 2 =tr(^(p)) = ^4|6 ifc | 2 , (5.42) 

j,k 

and simple àlgebra shows that 71 reverses T: 

K o F(p) = v 2 p . (5.43) 

Moreover, since only the restriction of 1Z to N is involved in the reversal and any oper- 
ation that reverses £ has the same restriction to N, we can conclude that any operation 
that reverses £ also reverses T. The converse is not true, however, because the canonical 
decomposition of T might map M to a set of orthogonal subspaces whose span is a proper 
subspace of N, in which case reversal of T would not entail reversal over all of N. 

The reversibility theorem proved in Sec. V A can be recast in another, very compact al- 
gebraic form. Suppose a quantum operation £ , with operator-sum decomposition consisting 
of operators Aj, can be reversed on a subspace M. Introduce the Hilbert-Schmidt inner 
product for operators, defined by 

(N, O) = tr(JVtO) . (5.44) 

This inner product allows us to define an adjoint of a superoperator, that is, any linear 
operator on operators. For example, the adjoint of £ is given by 

£\p) = Y J A) P A J . (5.45) 

3 

To see that this is an adjoint with respect to the Hilbert-Schmidt inner product, notice that 

(N,£(0)) = íJn^AjOA]] = tM^NA^O^ = {£\N),0) , (5.46) 

which is the required inner-product relation for an adjoint. The adjoint of an operation is 
generally not an operation, since it can be trace-increasing, but it is always a completely 
positive linear map. 

Let £m be the restriction of £ to the subspace M, as in Eq. (|4.5|) . Observing that 



£l o £ M ( P ) = J2 PmAÍAjPm p P M A)A k P M , (5.47) 



3,k 



we see that condition (|5.1|) is equivalent to the requirement that £* M o £ M be a positive 
múltiple of the identity operation on M. This requirement can be written as 
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Hi ° £m(p) =1 2 Pm pP M , (5.48) 
where the positive constant 7 2 is the trace of m 2 , that is, 

7 2 = EK fc r = E4 = ^EA?. (5-49) 

i,k j j 

Equation ( |5.48| ) is a necessary and sufíicient condition for reversibility of £ on M. We 
note that S M is generally not the reversal operator 71, for £j f has decomposition operators 
PmÀj = SdjPuU\ = ^JdjlljPj, whereas the reversal operation has decomposition operators 

shorn of the factor \J~dj. The omission is necessary so that 7Z is a trace-preserving operation. 

An equivalent way of writing condition ( |5.48j ) is to require that Em preserve, up to the 
constant 7 2 , the Hilbert-Schmidt inner product of operators defined on M; that is, 

(S M (N),S M (0)) = (N, £ M o S M (0)) = 7 2 (iV, P M OP M ) = -f 2 (P M NP M , P M OP M ) (5.50) 

for all operators N and O. 

VI. ERROR CORRECTION AND THE SECOND LAW OF THERMODYN AMICS 

Error correction — that is, reversal of an operation — decreases the entropy of a quantum 
system, so it is natural to inquire about the thermodynamic efficiency of this process. In 
this section we address the question of the entropy cost of error correction and show that 
error correction can be regarded as a sort of "refrigeration," wherein information about the 
system, obtained through measurement, is used to keep the system cool. Indeed, the method 
of operation of an error correction scheme is very similar to that of a "Maxwell demon," and 
the methods of analysis we use are based on those used to resolve that famous problem. As 
a prelude to our analysis, we review and extend the discussion of the Araki-Lieb inequality 
found in Sec. [III B| . 



A. Useful inequality 

The Araki-Lieb inequality |I2"-I4]| states that for two systems, 1 and 2, 



5(1) - 3(2) < 3 {12) . (6.1) 

To see this, introduce a third system, 3, which purifies 12. Subadditivity of the von Neumann 
entropy and the purity of 123 imply that 

S{1) = S(23) < S(2) + S{3) = S(2) + S{12) , (6.2) 



which gives the desired result. The inequality in Eq. ( |6.2| ) is the statement of subadditivity; 
equality holds if and only if systems 2 and 3 are in a product state, that is, 

p 23 = p 2 ® p 3 . (6.3) 
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By interchanging the roles of systems 1 and 2 in the above proof, the Araki-Lieb inequality 
can be written more generally as 

TO- 5(2)| < 5(12). (6.4) 

Suppose we apply the inequality ( |6.4| ) to a deterministic quantum operation V: 

\S(p R ')-S(pQ')\<S(p R V). (6.5) 

For a deterministic quantum operation, we have that S(p R ') = S(p R ) = S(p Q ); furthermore, 
it is always true that S(p RQ ) = S e (p,T>). Substituting these identities into the previous 
equation gives 

\S(p Q )-S(pV)\<S e (p,V). (6.6) 

A special case of this inequality is particularly useful in our entropic analysis of error cor- 
rection: 

S e (p, V) > S{ffl) - S(pQ') = -AS , (6.7) 

where AS = S(p® ) — S(jfi) is the change in the entropy of the system. From the equality 
condition (|6.3| ), we see that equality holds in the preceding equation if and only if 

p QE' =p Q'® p E> (6.8) 

These equality conditions are crucial to the following analysis of thermodynamically efhcient 
error correction. 



B. Reversal by a "Maxwell demon" 

Consider the error-correction "cycle" depicted in Fig. 1. The cycle can be decomposed 
into four stages: 

1. The system, starting in a state p, is subjected to a noisy quantum evolution that takes 
it to a state p n . We denote the change in entropy of the system during this stage by 
AS. In typical scenarios for error correction, we are interested in cases where AS* > 0, 
though this is not necessary. 

2. A "demon" performs a pure measurement, described by operators {-Bj}, on the state 
p n . The probability that the demon obtains result i is 

Vl = tv(B iP n Bt) , (6.9) 

and the state of the system conditioned on result i is 

Pi = B iP n B\/ Pi . (6.10) 

All error-correction schemes can be done in such a way that a measurement step of 
this type is included. 
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3. The demon "feeds back" the result i of the measurement as a unitary operation Vi 
that creates a final system state 



V iPi V} = ViBifTBM/pi 



(6.11) 



which is the same regardless of which measurement result was obtained. In the case 
of error correction this final state is the "corrected" state. 

4. The cycle is restarted. In order that this actually be a cycle and that it be a successful 
error correction, we must have p c = p. 

The second and third stages are the "error- correction" stages. The idea of error correction 
is to restore the original state of the system during these stages. In this section we show 
that the reduction in the system entropy during the error-correction stages comes at the 
expense of entropy production in the environment, which is at least as large as the entropy 
reduction. 



Noisy evolution 
(AS > 0) 



Restart cycle 
P c = P 



Measure {Bi} 
(H = B iP n Bj/ Pi 
Pi = tr(B lP n B\) 



"Feedback" result 
Pi p c = V^V} 



Environment 



"Demon'' 



Memory: i 



Figure 1. Error correction cycle. 



To investigate the balance between the entropy reduction of the system and entropy 
production in the environment, we adopt the "inside view" of the demon. After stage 3 
the only record of the measurement result i is the record in the demon's memory. To reset 
its memory for the next cycle, the demon must erase its record of the measurement result. 



Associated with this erasure is a thermodynamic cost, the Landauer erasure cost [19,20 



which corresponds to an entropy increase in the environment. The erasure cost of informa- 
tion is equivalent to the thermodynamic cost of entropy, when entropy and information are 



measured in the same units, conveniently chosen to be bits. Bennett 21] used the idea of 
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an erasure cost to resolve the paradox of Maxwell demons, and Zurek [22| and later Caves 



|23| showed that a correct entropic accounting from the "inside view" can be obtained by 



quantifying the amount of information in a measurement record by the algorithmic infor- 
mation content li of the record. Algorithmic information is the information content of the 
most compressed form of the record, quantified as the length of the shortest program that 
can be used to generate the record on a universal computer. We show here that the average 
thermodynamic cost of the demon's measurement record is at least as great as the entropy 
reduction achieved by error correction. 

In a particular error- correction cycle where the demon obtains measurement result i, the 
total thermodynamic cost of the error-correction stages is Jj + AS C , where 

AS C = S(p c ) - S(p n ) (6.12) 

is the change in the system entropy in the error-correction stages. What is of interest to us 
is the average thermodynamic cost, 

+ AS C ) = + AS C , (6.13) 

i i 

where the average is taken over the probabilities for the measurement results. To bound 
this average thermodynamic cost, we now proceed through a chain of three inequalities. 

The first inequality is a strict consequence of algorithmic information theory: the average 
algorithmic information of the measurement records is not less than the Shannon information 
for the probabilities Pi, that is, 

Em > W) = -T,pí 1 °£Pí • ( 6 - 14 ) 



Furthermore, Schack has shown that any universal computer can be modified to make 
a new universal computer that has programs for all the raw measurement records which 
are at most one bit longer than optimal code words for the measurement records. On such 
a modified universal computer, the average algorithmic information for the measurement 
records is within one bit of the Shannon information H. 

To obtain the second and third inequalities, notice that the corrected state p c can be 
written as 

P c = 5>*fal? = E ViBifTBÍV? = n(p n ) , (6.15) 

i i 

where 1Z is the deterministic reversal operation for the error-correction stages. The operators 
VíBí make up an operator-sum decomposition for the reversal operation. The probabilities 
Pi are the diagonal elements of the W matrix for this decomposition, 

Vl = ti{B lP n B\) = tT(V t B tP n B\VÏ) . (6.16) 

Thus we have our second inequality from Eq. ( |3.23| ), 

H(p)>S e (p n ,n) . (6.17) 

Equality holds here if and only if the operators VíBí are a canonical decomposition of TZ with 
respect to p n . We stress that different measurements and conditional unitàries at stages 2 
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and 3 lead to the same reversal operation, but they yield quite different amounts of Shannon 
information. 

The third inequality is obtained by applying the inequality (|6.7|) to 1Z and p n : 

S e {p n ,n) + AS c >0. (6.18) 
As Eq. (|5.37|) shows, equality holds here if the operators VíBí are a canonical decomposition 

of n. 

Stringing together the three inequalities, we see that the total entropy produced during 
the error-correction process is greater than or equal to zero: 

+ AS C > H(p) + AS C > S e (p n ,TZ) + AS C > . (6.19) 

i 

Stated another way, this result means that the total entropy change around the cycle is at 
least as great as the initial change in entropy AS, which is caused by the first stage of the 
dynamics. The error-correction stage can be regarded as a kind of refrigerator, similar to a 
Maxwell demon, achieving a reduction in system entropy at the expense of an increase in 
the entropy of the environment due to the erasure of the demon's measurement record. 

How then does this error-correction demon differ from an ordinary Maxwell demon? 
An obvious difference is that the error-correction demon doesn't extract the work that is 
available in the first step of the cycle as the system entropy increases under the noisy 
quantum evolution. A subtler, yet more important difference lies in the ways the two demons 
return the system to a Standard state, so that the whole process can be a cycle. For the 
error-correction demon, it is the error-correction steps that reset the system to a Standard 
state, which is then acted on by the noisy quantum evolution. For an ordinary Maxwell 
demon, the noisy quantum evolution restores the system to a Standard state, typically 
thermodynamic equilibrium, starting from different input states representing the different 
measurement out comes. 

Can this error correction be done in a thermodynamically efficient manner? Is there a 
strategy for error correction that achieves equality in the Second Law inequality (|6.19p ? The 
answer is yes, and we give such a strategy here. The proof of the Second Law inequality ( |6.19| ) 
uses three inequalities, J2iPdi > H, H > S e , and S e > —AS. To achieve thermodynamically 
efficient error correction, it is necessary and sufncient that the equality conditions in these 
three inequalities be achieved. 

We have already noted that Schack has shown that the first inequality, J2íPí^í > H{p), 
can be saturated to within one bit by using a universal computer that is designed to take 
advantage of optimal coding of the raw measurement records i. On such a universal computer 
the average amount of space needed to store the programs for the measurement records — 
that is, the encoded measurement records — is within one bit of the Shannon information 
H . Moreover, it is possible to reduce this one bit asymptotically to zero by the use of block 
coding and reversible computation. The demon stores the results of its measurements using 
an optimal code for a source with probabilities p^. Thus the demon stores an encoded list 
of measurement results. Immediately before performing a measurement, the demon decodes 
the list of measurement results using reversible computation. It performs the measurement, 
appends the result to its list, and then reéncodes the enlarged list using optimal block 
coding done by reversible computation. In the asymptotic limit of large blocks, the average 



28 



length of the compressed list of measurement results becomes arbitrarily close to H(p) per 
measurement result. 

The second inequality, H(p) > S e (p n ,TZ), can be saturated by letting the measurement 
operators Bi and conditional unitàries Vi be those defined be the those defined by the 
canonical decomposition of the reversal operation TZ. It should be noted that the optimal 
method of encoding the measurement records depends on the probabilities Pi, which in 
turn are ultimately determined by the initial state p. Thus the type of encoding needed to 
efficiently store the measurement record generally depends on the initial state p. For the 
canonical scheme for error correction, however, the probabilities Pi do not depend on the 
initial state p. 

The third inequality, S e (TZ, p n ) > — (S(p c ) — S(p n )), is satisfied by any error-correction 
procedure that corrects errors perfectly. Indeed, in Sec. |V B| we showed that the entropy 
exchange associated with any reversing operation is equal to the entropy reduction achieved 
by the reversing operation (see Eq. ( |5.37| )). An alternative demonstration that perfect 
error correction achieves the equality S e = —AS C begins by noting that at the end of the 
error-correction process RQ must be in a pure state-the initial state — and therefore the 
overall state must be a product p R ® <g> p EA (recali that E is the environment for the noise 
stage, while A is the ancilla for the reversal stage). Thus the condition p® A = p® £g> p A 
certainly holds. This is the equality condition (|6.8| ) for the Araki-Lieb inequality, applied 
to the reversal operation. Hence we have S e = —AS C for the reversal operation, and we 
conclude that any successful error-correction procedure automatically achieves equality in 
Eq. ( |6.18D . It would be interesting to see whether equality can be achieved in Eq. ( |6.18|) by 



error-correction schemes that do not correct errors perfectly. 



C. Discussion 



Zurek [25|, Milburn [26fl, and Lloyd f27j have analyzed examples of quantum Maxwell 



demons, though not in the context of error correction. Lloyd notes that "creation of new 
information" in a quantum measurement is an additional source of inefficiency in his scheme, 
which involves measuring <j z for a spin in a static 5-field applied along the z axis, in order 
to extract energy from it. If the spin is measured in the "wrong" basis — for example, if it 
is initially in a pure state not an eigenstate of a z — the measurement fails to extract all the 
available free energy of the spin, because of the disturbance to the system state induced by 
the measurement. In the case of error correction, something similar happens, but it is not 
disturbance to the system that is the source of inefficiency Instead, if the ancilla involved 
in the reversal decoheres in the "wrong" basis — that is, the measurement performed by the 
demon is not the one defined by the canonical decomposition of the reversal operation — then 
the Landauer erasure cost is greater than the eflicient minimum S e . This can be thought of 
a "creation of new information," due to "disturbance" of the ancilla, but the change in the 
system state is independent of the basis in which the ancilla decoheres. 

Error correction can be accomplished in ways other than that depicted in Fig. 1. The 
"inside view" of the preceding subsection, in which the demon makes a measurement de- 
scribed by some decomposition of the reversal operation, arises when the demon is decohered 
by an environment, the particular measurement being defined by the basis in which the en- 
vironment decoheres. If the demon is isolated from everything except the system and is 
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initially in a pure state, then its entropy gain is S e = —AS for the error-correction process. 
One can restart the error-correction cycle by discarding the demon and bringing up a new 
demon, the result being an increase in the environment's entropy by the demon's entropy S e . 
This way of performing error correction, which does not involve any measurement records, 
is equivalent to the "outside view" of the demon's operation. 

The "inside view" of the demon's operation, we stress again, arises if the demon's memory 
is "decohered" by interaction with an environment, the measurement record thus becoming 
"classical information." In this case the demon has the entropy H(p) of the measurement 
record, not just the entropy S e . Once this decoherence is taken into account, the different de- 
compositions of the reversal operation, corresponding to different measurements, constitute 
operationally different ways of reversing things, rather than just different interpretations of 
the same overall interaction. Keeping in mind the variety of decompositions of the reversal 
operation might lead one to consider a greater variety of experimental realizations, some of 
which may be easier to perform than others. As we emphasize above, a reversal in which 
the decohered measurement results correspond to a canonical decomposition of the reversal 
operation is the reversal method that is most efficient thermodynamically. 

VII. CONCLUSION 

In this paper we analyze reversible quantum operations, giving both a general informa- 
tion-theoretic characterization and a general algebraic characterization. Our results help in 
understanding quantum error correction, teleportation, and the reversal of measurements. 
By applying our two characterizations to a thermodynamic analysis of error correction, we 
show that the reduction in system entropy due to error correction is compensated by a corre- 
sponding increase in entropy of the rest of the world. Moreover, we show that error-correction 
schemes that correct errors perfectly can be done, in principle, in a thermodynamically ef- 
ficient manner. 
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